Quasi-Kondo Phenomenon due to Dynamical Jahn- Teller Effect 
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A mechanism of non-magnetic Kondo effect is proposed on the basis of a multiorbital Anderson 
model coupled with dynamical Jahn- Teller (JT) phonons. An electron system coupled dynamically 
with JT phonons has a vibronic ground state with double degeneracy due to clockwise and anti- 
clockwise rotational modes with entropy of log 2. When a temperature is lower than a characteristic 
energy to turn the rotational direction, the rotational degree of freedom is eventually suppressed and 
the corresponding entropy log 2 is released, leading to quasi-Kondo behavior. We discuss a simple 
relation between the "Kondo" temperature and the JT energy. 

PACS numbers: 75.20.Hr, 71.70.Ej, 75.40.Cx, 71.27.+a 



Forty years have passed since the pioneering work by 
Kondo to solve the problem of resistance minimum phe- 
nomenon Q . It has been revealed that the heart of his 
solution is the singlet formation from local magnetic mo- 
ment due to the coupling with conduction electrons Q, 
which is now called the Kondo effect Q . The Kondo- like 
phenomenon is expected to occur more widely, when a 
localized entity with internal degrees of freedom is em- 
bedded in a conduction electron system. Then, after the 
clarification of the Kondo effect with magnetic origin, lots 
of researchers have pursued a new mechanism of Kondo 
phenomena with non-magnetic origin. 

One possibility is to take into account explicitly or- 
bital degree of freedom of localized electron. First, Co- 
qblin and Schrieffer have derived exchange interactions 
from the multiorbital Anderson model Then, the 
concept of multi-channel Kondo effect has been devel- 
oped on the basis of such exchange interactions as 
a potential source of non-Fermi liquid phenomena ob- 
served in some /-electron compounds. Concerning the re- 
ality of the multi-channel Kondo effect, Cox has pointed 
out the existence of the exchange interaction in terms of 
quadrupole degree of freedom in a cubic uranium com- 
pound with non-Kramers doublet ground state p] . There 
have been continuous experimental efforts to detect such 
quadrupole Kondo effect. 

Another possibility for the non-magnetic Kondo be- 
havior is to consider local phonons. For instance, in 
an adiabatic approximation, Einstein phonons provide 
a double-well potential for an electron system, leading 
naturally to a two-level system, in which Kondo has first 
considered a possibility of non-magnetic Kondo behavior 
0. It has been shown that this two- level Kondo sys- 
tem exhibits the same behavior as the magnetic Kondo 
effect ||. Further efforts to include the low- lying levels 
of local phonon has been done intensively to unveil the 
quasi-Kondo behavior in electron-phonon systems 0. 

Recently, heavy-fermion phenomena in Pr- and Sm- 
based filled skutterudites have been experimentally found 
[TnL Hlj | . To understand such exotic properties of these 
compounds, the off-center motion of rare-earth atom in 



the pnictogen cage has been focused both from experi- 
mental |l2l| and theoretical viewpoints [tj • In particular, 
the off-center motion has been suggested to have degen- 
erate E g symmetry [12j. Inspired by this suggestion, we 
envision a scenario that Jahn- Teller (JT) effect plays a 
crucial role for the problem of Kondo-like behavior in 
filled skutterudites. From a conceptual viewpoint, this 
is a challenge to understand quasi-Kondo behavior by 
including both orbital and phonon degrees of freedom. 

In this Letter, we analyze a multiorbital Anderson 
model with the local coupling term between / electrons 
and JT phonons. Since such a coupling is believed to ex- 
ist in Pr-based filled skutterudites, we focus on a case for 
n=2, where n is the local /-electron number, correspond- 
ing to Pr 3+ ion. Using a numerical technique, we evaluate 
orbital susceptibility, entropy, specific heat, and phonon 
displacement. Then, we find that orbital fluctuations are 
enhanced in a vibronic state with double degeneracy due 
to rotational modes of clockwise and anti-clockwise di- 
rections, leading to residual entropy of log 2. When a 
temperature becomes lower than a characteristic energy 
to turn the rotational direction, the specific direction dis- 
appears due to the average over a long enough time and 
the entropy log 2 is released, leading to Kondo-like behav- 
ior. We also discuss a simple formula for the "Kondo" 
temperature of this phenomenon. 

Here we introduce the multiorbital Anderson model on 
the basis of a j- i co upling scheme for filled skutterudite 
compounds |13l Il5| . The main conduction band is 
constructed from p electrons of pnictogen 0] , expressed 
by a n with T7 symmetry in terms of the j-j coupling 
scheme. Note that / electrons in e u orbital with sym- 
metry is localized. Then, the model is given by 
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where is the dispersion of conduction electron, Cka 
is the annihilation operator of conduction electron with 
momentum k and spin a, and / 7<T is the annihilation op- 
erator of / electron on the impurity site with pseudospin 
a and orbital 7. The orbital index 7 distinguishes three 
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kinds of the Kramers doublets, two Tg ("a" and "b") 
and one Tr ("c"). V is the hybridization between con- 
duction and / electrons with T-j symmetry. Throughout 
this paper, we set V=0.05 and the energy unit is taken 
as half of the conduction bandwidth, which is about 1.4 
eV ^j- To adjust the local /-electron number as n=2, 
we appropriately set the chemical potential of / electron. 

The local /-electron term H\ oc consists of three parts as 
Hi oc =H C EF+H cc +H cph , where -ffcEF denotes the crys- 
talline electric field (CEF) potential term, H ee is the 
Coulomb interaction term among / electrons, and i/ op h 
indicates the coupling term between / electrons and lo- 
cal JT phonons. The CEF potential term is given by 
#CEF= J2-ya B if\afi<?, where B 1 is the energy level de- 
pending on orbital 7. Since the CEF term is already di- 
agonalized, it is convenient to introduce a level splitting 
A between T 7 and r 8 as A=i?r 8 — -E>r 7 - To reproduce 
the quasi-quartet situation with T% singlet ground and 
Ta triplet excited states for Pr-based filled skutterudites 
|18j . we prefer a small positive value for A as A=10 -5 . 

The Coulomb interaction term H cc is given by 
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where the Coulomb integral / in the j-j coupling scheme 
is expressed by Racah parameters, Eq, E\, and E2, which 
are set as Eq=5, Ei=2, and £^2=0.5. 

In the electron-phonon term -ff ph, the effect of E g rat- 
tling is considered to be included as the relative vibration 
of surrounding atoms. We remark that localized Tg or- 
bitals with e u symmetry have linear coupling with JT 
phonons with E g symmetry, since the symmetric repre- 
sentation of e u xe u includes E g . Thus, -ff e ph is given by 

H cph = g(Q 2 T x + Q 3 T z ) + (P 2 + P 2 )/2 

+ (uj 2 /2)(Q 2 2 + Ql) + b(Ql~2Q 2 2 Q 3 ), (3) 

where g is the clcctron-phonon coupling constant, Q2 and 
Q3 are normal coordinates for (x 2 — y 2 )- and (3z 2 — r 2 )- 
type JT phonons, respectively, P2 and P3 are corre- 
sponding canonical momenta, t x = Y,MLfb<r +fl ff f^), 
T z= J2a(fLf^ - flaha), w is the frequency of local 
JT phonons, and b indicates the cubic anharmonicity. 
Here the reduced mass for JT modes is set as unity. 
It is convenient to introduce non-dimensional electron- 
phonon coupling constant a and the anharmonic energy 
(3 as a=g 2 /(2uj 3 ) and /3=6/(2w) 3 / 2 , respectively. 

In this paper, to analyze the Anderson model, we ex- 
ploit a numerical renormalization group (NRG) method 
[20j, in which momentum space is logarithmically dis- 
cretizcd to include efficiently the conduction electrons 
near the Fermi energy. In actual calculations, we in- 
troduce a cut-off A for the logarithmic discretization of 
the conduction band. Due to the limitation of computer 
resources, we keep m low-energy states. Throughout this 



paper, we set A=5 and m=3000. Note that the temper- 
ature T is defined as T^A - ^ -1 ^ 2 in the NRG calcula- 
tion, where N is the number of the renormalization step. 
The phonon basis for each JT mode is truncated at a 
finite number -/V p h, which is set as _/V p h=25 in this paper. 

Since the dynamical JT phonons are coupled with or- 
bital fluctuations, we discuss orbital susceptibility to un- 
derstand the dynamical JT effect on electronic prop- 
erties. For the purpose, we evaluate the susceptibil- 
ity for corresponding multipole operator Ap 7 , where X 
denotes the multipole symbol and T is the irreducible 
representation with 7 to distinguish degenerate repre- 
sentations, due to the proper combination of three an- 
gular momentum operators J a (a=x, y, and z) for 
j=5/2. We note that orbital operators, t z and t x , corre- 
spond to 3g quadrupoles, 03 gu —(2J 2 — J 2 — Jy )/2 an d 
3gv =\/Z{J 2 — J 2 )/2, respectively, while another orbital 
operator, T y =-iJ2M Lfbcr - ft a M, is expressed by 
2u octupole, T2 U —Vl5JxJyJz/6, where the bar denotes 
the operation to take all possible permutations in terms 
of cartesian components [21|. Then, we evaluate three 
kinds of orbital susceptibility \ a for r a operators within 
the standard linear response theory. We also evaluated 
entropy S^mp and specific heat C- lmp for / electrons and 
average displacements, (Qi) and \J (Qf), for i—2 and 3. 

Before proceeding to the discussion on the effect of 
JT phonons, we briefly review the results for n—2 with- 
out JT phonons. As for details, readers should refer 
Ref. [l^. From the NRG results for multipole suscep- 
tibilities, we have found significant quadrupole fluctua- 
tions in addition to magnetic (dipole and octupole) ones 
at low temperatures. This result has been intuitively 
understood from the local electron configuration in the 
balance between Coulomb interaction and CEF potential 
[HI, Q [HI . For A>0, the local /-electron ground state 
is Ti singlet composed of two electrons. The first ex- 
cited state is Y± triplet formed by T7 and T% electrons, 
while the second excited state is T§ including a couple 
of Tg electrons. Even if Ti singlet is ground state, there 
exist T4 and T$ triplet states with small excitation en- 
ergy in Pr-based filled skutterudites. In fact, we observe 
residual entropy log 3 at low temperatures. Since these 
triplet states carry quadrupole moments, there remain 
quadrupole fluctuations in addition to magnetic ones. 

Now we include the effect of dynamical JT phonons, 
but without cubic anharmonicity. In Figs. 1(a) and 1(b), 
we show orbital susceptibilities, entropy, and specific heat 
for lu—0.3 and a=0.5. As mentioned above, due to the 
balance between Coulomb interaction and CEF poten- 
tial, triplet state exists near the Ti singlet ground state. 
Here we note that the T4 triplet state is JT active, since 
one electron is in the degenerate r 8 orbitals. Since the 
energy of the excited triplet state becomes lower due to 
the JT effect, we obtain a vibronic ground state with 
residual entropy log 2, in which Tg electron is tightly cou- 
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FIG. 1: (a) T\a and (b) Si mp and Ci mp vs. temperature for 
n=2, A=10 -5 , oj=0.3, a=0.5, and ,3=0. (c) Local electron 
potential in an adiabatic approximation for /3=0. 



pled with JT phonons. In the vibronic state for n=2, 
spin degree of freedom is suppressed in the non-Kramers 
states, while dynamical JT phonons enhance orbital fluc- 
tuations, as intuitively understood from the form of H cp h_. 
Then, three kinds of orbital fluctuations for t x , t v , and 
r 2 significantly remain in the low-temperature region. 

Let us here discuss the nature of the vibronic state. For 
the purpose, it is useful to see the electron potential in an 
adiabatic approximation, as shown in Fig. 1(c), although 
in actuality, the potential is not static, but it dynami- 
cally changes to follow the electron motion. For ,9=0, 
the potential is continuously degenerate along the circle 
of the bottom of "sombrero" . Thus, there should exist 
double degeneracy for rotational modes along clockwise 
and anti-clockwise directions, as we observe the entropy 
of log 2 in Fig. 1(b). When a temperature becomes lower 
than a characteristic energy T* , which is related to a 
time scale to turn the direction of rotational JT modes, 
the residual entropy log 2 should be eventually released, 
leading to Kondo-like behavior, since the specific rota- 
tional direction disappears due to the average over the 
long enough time. However, we cannot see such behavior 
in the present temperature range, because T* is consid- 
ered to be very low. 

In order to detect the quasi-Kondo behavior, we in- 
crease T* by decreasing the excitation energy concerning 
the rotational JT mode. Such an energy is closely related 
to the static JT energy Ejt in the adiabatic approxi- 
mation. In order to decrease effectively Ejt by keeping 
the coupling strength, we include cubic anharmonicity. 
As shown in Fig. 2(a), for (3<0, three potential min- 
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FIG. 2: (a) Local electron potential in an adiabatic approxi- 
mation for /3<0. (b) T\a (a—x, y and z), (c) Simp and Ci mp , 
and (d) {Qi) and y / {Qf) (i=2 and 3) vs. temperature for 
n=2, A=1CT 5 , lu=0.3, a=0.5, and /3=-0.006. (e) Schematic 
views for displacements for T>T* and T<T*. 



ima appear in the bottom of the potential. Since the 
rotational mode should be due to the quantum tunneling 
among three potential minima at low temperatures, the 
frequency is effectively reduced in the factor of e ~ SE / u ^ 
where SE is the potential barrier. Then, we expect to 
observe the quasi-Kondo behavior even in the present 
temperature range, when we include the effect of cubic 
anharmonicity. 

In Fig. 2(b), we show orbital susceptibilities for w=0.3, 
a=0.5, and (3=— 0.006. In the present unit, SE is about 
100 K. For 10 -7 ;$T<0.01, we observe orbital properties 
similar to those in Fig. 1(a), since the vibronic ground 
state also appears in this temperature range. However, 
with decreasing temperature, we observe the lift of de- 
generacy in orbital fluctuations and the gradual decrease 
of T\x and T\ y around at T=10~ 8 . Then, only the 
t z fluctuation remains at low temperatures. As shown 
in Fig. 2(c), at the crossover temperature T*, we actu- 
ally observe the release of entropy log 2 and a clear peak 
in the specific heat, suggesting the Kondo-like behavior, 
even though T* is still low. 
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FIG. 3: Characteristic temperature T* vs. Ejt- Dashed line 
is a guide for eyes. Size of the symbol indicates the order of 
the precision of T* in the calculation. 



tiorbital Anderson model coupled with dynamical JT 
phonons. We have observed the remarkable Kondo-like 
phenomenon induced by dynamical JT effect and dis- 
cussed a simple formula for "Kondo" temperature. As for 
reality in actual materials [24^. rattling motion in filled 
skutterudites provides dynamical JT phonons coupled to 
degenerate orbitals. Thus, it is interesting to pursue pos- 
sible relevance of the quasi-Kondo phenomenon to exotic 
heavy fermion behavior in filled skutterudites. 
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In order to understand what happens in this quasi- 
Kondo phenomenon, let us show the temperature de- 
pendence of average displacements. For T*<T^0.01, we 
observe y/{q^=^{Qf)^0 and (g 2 > = (Q 3 >=0. Namely, 
there is no average displacement and the vibration is al- 
ways isotropic, as schematically shown in Fig. 2(d). How- 
ever, after the release of log 2 entropy for T<T*, we ob- 



=|(Q 3 )M0 and v/(Q|)=(Q 2 )=0, indicating 



tain yj (Qf) 

that the vibration becomes anisotropic due to the selec- 
tion of the Q3-type mode, as shown in Fig. 2(d). Then, 
the degeneracy in orbital fluctuations is lifted and only 
the t z fluctuation remains at low temperatures. 

Finally, let us intuitively discuss a simple formula for 
the characteristic temperature T*. For such a purpose, it 
is convenient to consider the exchange interaction J for 
this case. Since J should be related to the process of dis- 
appearance of JT distortion, it is roughly expressed as 
JocV 2 /Ejt, where Ej T is given by Ej T =au>+&\l3\a 3 / 2 . 
Here we recall the expression for the Kondo tempera- 
ture, 2k=Wexp(— 1/Jpo), where W is the bandwidth of 
conduction electrons and po is the density of states at 
the Fermi level. Then, we deduce the relation for T* 
as logT*oc— Ejt- In order to test the validity of this 
relation, in Fig. 3, we plot logT* vs. Ejt for several 
parameter sets, where T* is defined as a temperature 
showing a peak in Ci mp . Within the errorbars, we can 
conclude that the relation of logT* cx — -Ejt holds for 
the quasi-Kondo phenomenon. A similar expression for 
the Kondo temperature was obtained in the context of 
odbital Kondo effect jl^. Note, however, that this re- 
lation cannot be applied to all parameter regions, since 
it is simply deduced in the adiabatic approximation for 
JT phonons. We also note that the effect of Coulomb 
interactions is ignored, but it appears as an additional 
effect for the renormalized JT energy j^, and thus, 
the expression for T* does not change qualitatively. 

In summary, we have numerically analyzed the mul- 
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